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We consider the touchdown of an elastic sheet as it adheres to a wall, which has a dynamics that 
is limited by the viscous resistance provided by the squeeze flow of the intervening liquid trapped 
between the two solid surfaces. The dynamics of the sheet is described mathematically by elastohy¬ 
drodynamic lubrication theory, coupling the elastic deformation of the sheet, the microscopic van 
der Waals adhesion and the viscous thin film flow. We use a combination of numerical simulations of 
the governing partial differential equation and a scaling analysis to describe the self-similar solution 
of the touchdown of the sheet as it approaches the wall. An analysis of the equation satisfied by the 
similarity variables in the vicinity of the touchdown event shows that an entire sequence of solutions 
are allowed. However, a comparison of these shows that only the fundamental similarity solution 
is observed in the time-dependent numerical simulations, consistent with the fact that it alone is 
stable. Our analysis generalizes similar approaches for rupture in capillary thin him hydrodynamics 
and suggests experimentally verihable predictions for a new class of singular hows linking elasticity, 
hydrodynamics and adhesion. 


The mathematical description of huid hov^^ often leads to diverging velocity gradients in bulk hows or 
geometric singularities in free surface how^^JEl^ with implications for many applications such as microhuidics, 
coating, food processing etc. Here, we describe the elastohydrodynamic thin him how associated with van 
der Waals-driven touch-down between an elastic sheet and a solid substrate, whic h is i nspired by the rupture 
of capillary thin huid hlms on substrate^^®^ that reveal a hnite-time sin gularitjl^liSI 

Our analysis focuses on small scale hows wherein inertial ehects from both the huid and the sheet may 
be neglected. Furthermore, we consider the touchdown of a long and wide elastic sheet with a length L, 
which is initially separated from a solid substrate by a thin viscous huid him of height h{x,t). The sheet 
is destabilized and attracted to the substrate by a van der Waals adhesion potential (Fig. la). Mechanical 
equilibrium requires that it the shape of the sheet is determined by the equation for the pressur^i^p(i, t) 

A 

p{x,t) = Bhi±s,x{x,i) - . , (1) 

3h'^[x, t) 


where Bhxxxx{x,i) is the elastic bending pressure due to long wavelength deformations of the sheet with 
bending stihness B [N-m], and - . is the van der Waal adhesion pressur^i^ with A [N-m] the Hamaker 

constant. Here, and elsewhere, the subscript denotes derivative i.e. {■)x = d{-)/dx. Furthermore, we have 
assumed that the fluid film is very thin, so that its thickness h{x,t)/L <C 1 and that it has a small slope 
hx{x,t) <C 1, and finally that there are no inertial effects associated with the motion of the elastic sheet. 

As the sheet approaches the substrate the intervening fluid is squeezed out, leading to non-uniform de¬ 
formations of the sheet. In the small gap approximation, viscous forces dominate over inertial effects, so 
that we may use to the lubrication approximation for the description of the fluid motiorP^ that yields an 
evolution equation for the film height h{x, given by 
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p is the fluid viscosity and we have assumed that there is no-slip of the fluid at both surfaces. Using Q, 
we find that the resulting sixth order nonlinear partial differential equation (H) couples the dynamics of the 
elastic sheet, the viscous forces in the lubricating film and the destabilizing intermolecular adhesion pressure. 
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We scale the variables appearing in (|^ as; a; = t'x = xho{B/A) i , h{x, t) = hoh{x, t) and i = rt = 

where r = m jg viscous time scale, and introducing these scaled variables into ([^ gives us the 

parameter free dimensionless elastohydrodynamic lubrication equation, 

ht{x,t) = (h^ix,t)hxxxxxix,t)- ■ (3) 

To complete the formulation of this initial value problem, we complement ^ by an initial condition 
ho = h{x,t = 0) = 1 — 0.05 X (1 + cos(a;)) corresponding to a film with a constant height that has a small 
blip at a; = 0, and six boundary conditions at the ends. Since the dynamics of interest in the touchdown 
region is far from the boundary, we use symmetry conditions at the two ends of the sheet, resulting in 
vanishing derivative^ i.e. hx = hxx = hxxx = 0. Note that these boundary conditions does not affect the 
dynamics near the touchdown at x = xc, t ^ tc- 

Before solving numerically it is of interest to briefly consider the dispersion relation of the linearized 
version of ([^. Substituting h{x,t) = 1 + eexp(At + ikx) with e ^ 1 in (|^ with k the wave number, A the 
growth rate, we find that A = fc^(l — /c^). Thus, we see that the elastic thin film is only unstable for long 
wavelength perturbations fc < 1, similar to what is known in capillary dominated flows. This suggests that 
when we solve © numerically a minimal domain size x > 2TTjkc is required, with kc = 1 the critical wave 
number. 

We solve © using a second order finite difference discretization in space and a Gear methocP^ for the 
time stepping. In Fig. la we plot the evolution of the shape of the sheet and see that the solution appears 
to adopt a self-similar shape in the neighborhood of the touchdown point x = xc = 0,t = tc, where the 
contact time tc is defined as the last time point before h < 0, just before the numerics breaks down. 

Inspired by the self-similar shape of the sheet as —>■ 0 in the neighborhood of x = xq, t = tq, we look 
for a self-similar solution of © to understand the scaling behavior of the sheet in the immediate vicinity of 
touchdown. Adopting the similarity ansatz 

Hx, t) = {tc - tHH{r]), r] = (4) 

with H{r]) the scaled height, rj the similarity variable, and inserting Q into ([^ gives 

{t - tcT~^{aH{r]) - /3f]H'{r])) = (^{t - tc)‘^°‘~^I^H{r]fH'""{T]) - {t - (5) 

with H'{r]) = dH{ri) /drj. Equating the exponents appearing on both sides of the equation, we find that 
a — 1 = 4a — 6/3, a — 1 = —2/3 so that a = /3 = 1/3 and are the same in axial symmetrjli^. Inserting this 
result in ([^ we find that the similarity variable H(rj) satisfies the ordinary differential equation 

i (77(77) - vH'iv)) = ■ ( 6 ) 

To test the postulated self-similar scaling ansatz Q with a = /3 = 1/3, we re-plot the results obtained 
by solving the partial differential equation ([^ in Fig. la by rescaling h{x,t) using the ansatz above. In 
Fig. lb, we see that the solution of ([^ converges rapidly to a self-similar form H{r]) over many orders of 
magnitude in 77 , confirming the validity of our similarity hypothesis. 

The collapse of the time-dependent numerical results from ([^ by using the similarity transform in Q 
shows the self-similar shape of the sheet. Since this solution is also described by ©. a natural question is 
if we can solve this equation directly and determine the self-similar shape associated with the touchdown 
event. This requires that we supplement ([^ with a set of boundary conditions. Noting that the numerical 
solution of ^ in Fig. la is symmetric around the contact point 77 = 0 suggests that odd derivatives of the 
solution vanish, i.e. H'{0) = H”'{Q) = 77"'"(0) = 0. Furthermore, since the sheet height chan ges rapidly 
near the contact point i.e. x = xc, t —)■ tc, the far-held is quasi-static by comparison i.e. ht « (pESl. Thus, 
we expect that the far-held solution is independent of the evolution of the singularity as t —> tc- This 
translates into a far-held Robin boundary condition that in similarity variables reads as 


i (77(77) - 7777'(77)) = 0, I 77 1^00. 


( 7 ) 
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(b) (c) 

FIG. 1 . (a) Numerical simulations of ([^ for the dynamics of tonchdown showing the evolntion of the sheet at five 
different time points t = [0, 6.50, 7.98, 8.08, 8.09]. The sheet is initialized with ho — h{x, 0) = 1 — 0.05 x (1 + cos(a:)), 
and is seen to adopt a linear profile as t ^ tc for x £ [—2.5 2.5). Our simnlations allow us to measure the contact 
point xc = 0 and contact time tc = 8.0899. (b) The shapes obtained by solving § may be collapsed onto a universal 
self-similar shape by using the rescaled variables H{ri) — h{x,t)/{tc — t)^ ,r] = {x — xc)/{tc — , over more than 

two orders of magnitude in H{ri) and rj. The inset in the lower right shows the self-similar collapse for r} £ [—20 20]. 
(c) The rescaled height H{rj) obtained by solving the PDF (§ and the similarity ODE ([^ 7 ) (with m = 1 ) shows 
that the height profiles adopt a shape H{rj) = a\r] as ] rj ]—^ oo (o). The shaded area »? < 1 illustrates the inner 
region near the contact point where the shape of the sheet arises from the balance between viscous forces, van der 
Waals adhesion pressure and elastic bending resistance. 


We immediately see that the far-held solution is satished by H{r]) = aij + C as \ ry |—)■ oo, and the integration 
constant C = 0 for consistency. Plotting the results of the numerical solution of § shown in Fig. la in terms 
of the rescaled height H{ri) as t ^ tc in logarithmic coordinates (Fig. Ic) conhrms that H{t]) ^ rj over a wide 
range of 77 , so that the constant a is then interpreted as an asymptotic matching condition for the far-held 
associated with ([^, requiring us to prescribe an additional boundary condition. Since H{r]) ^ arj as rj ^ 00, 
this leads to the following boundary conditions at inhnity; H'{oo) = a, H”{oo) = H"'{oo) = H""{oo) = 0, 
which together with the symmetry conditions at 77 = 0 complete the formulation of the boundary value 
problem for H{rii) that satishes ( 6 ). 

We solve the self-similar ODE (6-7) with the hnite-difference-based boundary value solver BVP4C in 
Matlat(i^ on a hnite size domain ry = 50, with the boundary conditions; H'{0) = H'''{0) = = 

H"{50) = i7"'(50) = H'"'{50) = 0, 77 ' (50) = a, noting that our results are insensitive to changes in the 
domain size. We hnd that we get a discrete family of solutions for ([^7) similar to those seen in the problem 
of capillary ruptur^. The subscript m is used to distinguish the different discrete solutions of 77(ry)and a 
with 


77(77) — 77771(^)7 a — CLm ' 


( 8 ) 


The numerical procedure starts with an initial guess for the height at contact 77o(0), together with a guess of 
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m 


Hm{0) H"(0) Rr(0) 

PDE (3) 

- 

0.708 

0.324 

-0.183 

1 

0.625 

0.719 

0.321 

-0.186 

2 

0.290 

0.452 

0.031 

0.190 

3 

0.193 

0.338 

0.051 

-0.099 

4 

0.147 0.283 

0.0215 

0.068 

5 

0.119 

0.245 

0.0223 

-0.034 


TABLE I. Discrete solutions of the self-similar ODE (6-7) showing the values am, -ffm(O), Hm(Q), H'm {Q). Comparison 
of these values with those obtained from the solution of the PDE § show that it is consistent only with the 
fundamental self-similar solution Hi{ri). 


the initial solution Ho{r]) = Ho{0) + aor]. We systematically vary the initial guesses Ho{0) and oq, allowing 
us to identify the first thirteen solutions m = 1, 2,13 for Hmirf) satisfying §7), see Fig. 2a . 

Comparison of these solutions of ODE ([^7) with the solution of the time-dependent partial differential 
equation ([^ shows that only the fundamental solution associated with to = 1 agrees with the PDE in ([^, 
as shown in Fig. 2a. 

To understand why only the fundamental solution to = 1 is consistent with the time-dependent solution of 
^ we qualitatively evaluate the stability of these self-similar solutions, using Hmiv) as the initial condition 
for the time dependent PDE In Fig. 2b, we show the result of our simulation of ^ with h{x, 0) = Hmij]), 
TO = 3. The time and position of contact shifts and similar results arise for all values of to > 1. Similar 
observations have been made for capillary thin film ruptur^. Rescaling the results with (4) shows that the 
solution converges to that associated with the shape of the fundamental solutions to = 1 in the vicinity of 
touchdown (Fig. 2c). While a detailed stability analysis of the problem is necessary to classify the different 
unstable modes, we postpone this to the future and instead give a simple heuristic explanation next. 

Starting with the different discrete solutions of the similarity equation Hmiff), we see that for to > 1 
there are small pressure oscillations that are defined by the bending pressure Hmiv) in the neighborhood of 
77 = 0. As TO increases am decreases and there is a subsequent flattening of the film (Table 1, Fig. 2a). The 
number of maxima/minima increase linearly with to as can be seen in Fig. 2d. Since these correspond to 
short wavelength modes of the sheet, they are damped out much more strongly than the longest wavelength 
mode corresponding to Hi{r]) which is balanced by the van der Waals adhesion pressure H~^{r]) and the 
viscous stress, and this self-similar solution is the only dynamically stable one that survives. 

We conclude with a few remarks on the approximations inherent in our approach. As the minimal 
height becomes comparable to the molecular dimension « 2.5A, the continuum description that we have 
used breaks down. Our simulations show that the self-similar shape is adopted already at h{x,t) « 0.3 
(Fig. lb) and choosing a typical value for the initial height Hq = 300nm gives a height « lOOnm which is far 
from this molecular height and within experimental realm using total internal reflection microscopy imaging. 
During touchdown, the minimal height hmin(t) and the horizontal deformation length £{t) = hminit){B/A)* 
becomes dynamic. If £{t) is comparable to the thickness of the sheet b ~ hmin{t){B/A)^, three-dimensional 
elastic effects may become relevant. However, since our small slopes approximation holds throughout time 
as hx hmin{t)/£i£) ^ 1 with B ^ A e.g. A = 10“^°N-m, B = IN-m the plate theory that 

we have used here remain reasonable. Finally, as the sheet and fluid velocity diverges near the singularity, 

this can lead to inertial effects, particularly as the reduced Reynolds number Re = ~ {tc — t)~^ 

is also singular as t —)■ tc- By assuming Re « 0(1), we get a scaling estimate for the height at which 


inertia becomes significant hj 


, 3 
pA2 


e.g. 


77 = 10 ^ Pa-s, p = 10^ kg/m^, giving hp < 2.5A and the 


12B2fi 

approximation of neglecting inertial effects remains valid. 

Our analysis reveal a new class of singular flows linking elasticity, hydrodynamics and adhesion, relevant 
to contact between a thin elastic sheet and an adherent surface. A combination of numerical simulations and 
similarity analysis shows that the dynamic height and horizontal deformation of the sheet in the process of 
contact show a power-law form. Analysis of the governing equations in similarity coordinates further allows 
us to determine the universal self-similar shape of the elastic sheet in the neighborhood of touchdown. A 
qualitative analysis of the discrete set of self-similar film shapes that satisfy the similarity ODE shows that 
the solution with the least oscillatory bending pressure is the one that is consistent with the time-dependent 
numerical simulation, while others with short-wavelength pressure oscillations are strongly damped. A 
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FIG. 2. (a) Comparison between the self-similar shape of the sheet predicted by the time-dependent PDF © and 
the first thirteen solutions Hm{ri) \ m £ [1 — 13] of the similarity ODE ©7). We note that only the first fundamental 
solution H\{r]) of (6-7) is consistent with that of the PDE © (o). (bjNumerical simulations of the PDE © using 
the initial condition associated with the self-similar shape ho = h{x,t = 0) = 773(h) shows the evolution of the 
shape of the sheet for four different time points t = [0, 0.28, 0.336, 0.341]. We find that the numerically measured 
contact time is tc = 0.3412 and the contact points are xc = [—1.48,1.48]. (c) By using the rescaled variables 
77(h) — h{x,t)/{tc — t)3 ^r] = {x — xc)l{tc — 7)® with xc = 1.48 corresponding to the contact on the right side, we 
see that the shapes of the sheet obtained in Fig. 2a collapse onto the universal shape of the fundamental solution 
of (6-7) 77i(h) (shown in red). The deviations away from the contact point are due to the effects of the second 
touchdown at Xc ~ —1.48. (d) To distinguish the different discrete solutions Hmijj) satisfying (6-7) we plot the 
scaled bending pressure 77((,"(h), and see that the number of maxima/minima scales with m. Thus, for m > 1, the 
wavelength of the pressure oscillations decreases as m increases; these solutions are therefore strongly damped out 
in comparison with the fundamental solution of H\(rf) (6-7), which alone survives as we approach touchdown. 


natural next step is to experimentally test the regime of applicability of our results. 
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